Green function approach for scattering quantum walks 



O 
(N 

o 

Q 



Oh- 



> 



X 



F. M. Andrade^'El and M. G. E. da Luz^'Q 
^ Departamento de Matemdtica e Estatistica, Universidade Estadual de Ponta Grossa, 84030-900 Ponta Grossa-PR, Brazil 
^ Departamento de Fisica, Universidade Federal do Parana, CP. 19044, 81531-980 Curitiha-PR, Brazil 

(Dated: December 13, 2011) 

In this work a Green function approach for scattering quantum walks is developed. The exact 
formula has the form of a sum over paths and always can be cast into a closed analytic expression 
for arbitrary topologies and position dependent quantum amplitudes. By introducing the step 
and path operators, it is shown how to extract any information about the system from the Green 
function. The method relevant features are demonstrated by discussing in details an example, a 
general diamond-shaped graph. 
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I. INTRODUCTION 

Generally speaking, quantum walks (QW) represent 
unitary evolutions taking place in discrete spaces - 
graphs - for which typical basis states are localized. 
There are several ways to formulate QW, either consid- 
ering time as a continuous (CTQW) 0] or a discrete 
variable. In the latter case, the two major formulations 
are the (a) coined QW (CQW), based on inner "coin" 
states (see, e.g., Q), and (b) scattering QW (SOW), re- 
lying on the idea of multi-port interferometers 0)1^1 ■ The 
continuous time and coined QW are directly related as a 
limit process 0, whereas the CQW and SQW have been 
shown to be unitarily equivalent in arbitrary topologies 

Quantum walks originally emerged 8] from the interest 
to construct and understand quantum analogs of classi- 
cal random walks (CW). But soon it was realized they 
also would constitute powerful tools in quantum compu- 
tation , specially given that QW can represent univer- 
sal quantum computation primitives [lo| . In fact, for a 
long time CW have been used to solve different computa- 
tional problems [ll|. Thus, the connections between the 
quantum and classical walks I12h14| |. allied to the par- 
ticular features of the former |l5| , actually point to the 
potential usefulness of QW in building algorithms which 
are much faster and robuster [l^, than their classical 
counterparts. As r epr esentative examples we can cite the 
Grover algorithm isj (simulated through QW lo*]) for 
searching of unsorted database, the element distinctness 
algorithm [23| , the detection of marked elements [IJ , the 
computation of orders of solvable groups [l^, and the 
quantum Fourier transform [23| . Moreover, even prob- 
lems like the energy transport in biological systems can 
be analyzed by means of Q W [13] . 

A key aspect in such class of systems is the quantum 
interference between the possible "paths" (see next Sec- 
tion) along the evolution (25i - i28i] . It leads to a dynam- 



ics that generally spreads much faster than CW fis'l (al- 
though in certain situations anomalous sub-diffusive be- 
havior may also emerge As a consequence, one get s 
exponentially faster hitting times from QW P, HO, ISll . 
one of the reasons why QW are particularly suitable [32| 
to solve searching problems [27j, [S^. Also, different dif- 
fusion processes, from ballistic to Anderson localization 
[33, HI] , are possible when decoherence is included. 
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Since interference is fundamental to explain different 
phenomena observed in QW [HI (including many of the 
applications mentioned above) it is desirable to have a de- 
scription emphasizing the path-like character of QW. In 
this respect Green function methods are particularly use- 
ful (37H39I . Then, here we develop a full Green function 
approach for QW in arbitrary topology and for position 
dependent quantum amplitudes. For so, we assume the 
very appropriate discrete scattering formulation, SQW. 
We should observe there are few interesting works (e.g., 
Refs. (40l - |42j ) addressing the classification of trajectories 
in QW. They, nevertheless, are based mostly on combi- 
natorial analysis to compute all the possible final states 
at a time t = n. Our proposed construction thus is much 
closer to the idea of "history" of trajectories in the Feyn- 
man sense j43i] . 

The paper is organized as the following. In Sec. II 
we review the scattering formulation for quantum walks, 
also making few useful parallels with the classical case. 
By direct mapping one-dimensional QW to a related type 
of problem, ID point interaction lattices, in Sec. Ill we 
are able to write the exact Green function G in the form 
of a sum over paths. Moreover, we discuss how such for- 
mula can be summed as a closed analytical expression. 
In Sec. IV the ID construction is extended to complete 
arbitrary topologies. By defining the step and path op- 
erators, we show in Sec. V how to extract any system 
relevant information from the exact expression for G. In 
Sec. VI we illustrate the features of the present approach 
analyzing in details a particular example, the diamond- 
shaped graph. Finally, we present the conclusion in Sec. 
VII. 
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FIG. 1. The "Hilbert lattice" associate to SQW in ID. For 
each site it is defined appropriate scattering quantum ampli- 
tudes (here illustrating the phase convention in [J). 



II. A BRIEF REVIEW ON THE SCATTERING 
FORMULATION FOR QUANTUM WALKS 

To better understand the main ideas underlying the 
definition of quantum walk models, and thus to develop a 
Green function approach, here we review QW scattering 
formulation [i] on the line (ID). The case of more general 
topologies wiU be discussed in the next Sections. 

So, consider a helpful framework for SQW: to view 
their evolution as a dynamics defined on a ID "Hilbert 
lattice", depicted in Fig. 1. Notice, however, it does not 
necessarily represent a spatial structure since the states 
(assumed on the bonds) do not need to be position eigen- 
vectors. Under this picture, the lattice characteristic pa- 
rameter \a L ^ Aj — 1, just the spacing between two 
consecutive sites of the Hilbert lattice. Along each bond, 
joining the sites j and j + 1 (Fig. 1), we have two pos- 
sible states, I -I- l,j -I- 1) and | — l,i). Then, each ba- 
sis element, |cr, j), is labeled by two quantum numbers. 
The first, a, sets the "direction" (±1) along the lattice. 
We mention that although fully equivalent, the present 
is slight different than the common SQW construction in 
the literature. 

The discrete time evolution is given by the one step 
unitary operator U, such that |^'(n + 1)) U\'i'{n)). 
For U, we consider the translation (St) and reversion- 
translation (Srt) operators 

SrlcrJ) = + cr), St^ctJ) = \a,j - a), 
SrtWJ) = Srt^ctJ) = \ - (J,j - <t), (1) 

with both unitary and Srt^ ~ 1. We also define T and 
I 



R, for which the basis states \(J,j) are eigenstates, or 

T\<tJ) ^t„j\a,j), R\a,j) = r„J(T,j). (2) 
If we impose now 

+ ij i-hj = + rlj^j t+ij = 0, 



then, the unitary time evolution reads (0 < 7 < 27r) 

Uij) = eMn] {St T + Srt R) . (4) 

The term exp[z7] is associated to the translation between 
neighbor sites (Aj = L = 1), relevant to proper describe 
stationary scattering solutions 44 1. 

Provided Eq. holds, there is a freedom to choose 
the coefficients Vj and tj. For instance, by setting (0 < 
Pj < 1 and < (fij, ipj < 2tt for any j) 



t„ 



<jJl - p'^.exp[iaipj], (5) 



one gets = — r^i j and t-ij = t*^^ j, just the con- 

vention used in Q (Fig. 1). 

The dynamics in Eqs. (Il])-(l4]) in fact represents an ex- 
tended quantum version of a more simple classical ran- 
dom walk. Each time the classical walk needs to choose 
a new direction, it uses the same probabilities (P and 
1 — P) to decide between right and left. By allowing in 
Eq. (O p and the phases to depend on j, we arc implic- 
itly assuming position dependent distribution functions 
for the direction choices. Obviously, by setting a same p, 
if and (f) for any j we recover the usual case. 

Finally, as it stands the above model is deterministic in 
the quantum mechanical sense: any initial state |^'(0)), 
after n time steps, is uniquely determined by the always 
weU defined state ?7"|^'(0)). Thus, stochasticity (i.e., 
classical randomness) can enter into the problem only 
through measurements, when we determine the system 
location along the Hilbert lattice. In fact, 



P.j(n) = |(j,a|vl/(n))|^ 



(6) 



is the probability to be in the quantum state (or in the 
present lattice language "position and direction") j,a at 
time n. So, projection is an essential ingredient in QW. 

As a simple example, consider the initial state |5'(0)) = 
I -I- 1,0). Under U one has after n = 3 time steps 



|*(3)) = exp[3i7] 1^+1,0^+1,1^+1,2 1 + 1,3) + (r+i,or_i,_it+i,o + t+i,o?'+i,i'^-i,o)l + 1,1)+ r+i,ot_i,_ir_i,_2 | + 1,-1) 
+r+i,o<-i,-it-i,-2| - 1, -3) + (7-+i,or-i,-ir+i,o + i+i.o'"+i,it-i,o)| - 1, -1) + <+i.o<+i,i'"+i:2| - 1, 1)}- (7) 



Thus, the system probability to be found, say, in | + 1, 3) there is only one possible "path" ending up in | + 1, 3). 
is |i+i,oi+i,it+i,2p- Note that for three time steps. Hence, the modulus square of the quantum amplitude 
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FIG. 2. Schematic association between QW Hilbert (space) 
lattice and a generalized Kronig- Penney configuration lattice. 
Each site j corresponds to a point interaction ai x — jL (of 



reflection and transmission amplitudes rj^' and 



associated to such path yields the sough probability. 
On the other hand, there are two possible paths lead- 
ing to I — 1,-1). They correspond to the amplitudes 
r+i,o''-i,-i''+i,o and i+i,or+i,ii_i,o (cf. Eq. ©). But 
contrary to CW, where the total probability is the sum 
of the individual probabilities of each trajectory, here the 
quantum interference character of the walk demands that 
P_i,_i(3) = |r+i,o?'-i,-i?'+i,o +t+i,or+i,it^ifi\'^. 



III. A GREEN FUNCTION APPROACH FOR 
QUANTUM RANDOM WALKS 

Here we develop a Green function approach for the 
SQW in the previous Section, proceeding in three steps, 
(a) First, we construct a mapping from our QW to a 
ID generalized Kronig- Penney lattice [1^, for which we 
can calculate the exact energy- dependent Green function 
G. (b) Then, we discuss which are the mapped system 
appropriate configurations in order to match the original 
problem, (c) Finally, we show how the obtained G gives 
the quantum walk sought dynamics. We leave to the 
next Section the extension of the ID results to the case of 
more general topologies, namely, QW on arbitrary graph 
structures. 



A. The mapping 

As already emphasized, the quantum walk does not 
necessarily represent any dynamics on a concrete physical 
lattice. Nevertheless, for our purposes it is very useful to 
associate the quantum walk Hilbert (space) lattice - and 
its underlying "kinematics" [i^ - to that of an usual 
continuous ID quantum scattering problem. 

In Fig. 2 we show schematically the correspondence 
between the model of Fig. 1 with a generalized Kronig- 
Penney lattice of equally spaced arbitrary point interac- 
tions [381 1 ■ i.e., zero range potentials which extend the 



usual delta function [4g]. Each point interaction (at 
X = ±jL, j = 0, 1, . . .) is entirely characterized by the 
quantum amplitudes r^^\k) and t^j^\k). The super- 
script -|- (— ) stands for the reflection or transmission 
of a plane wave of wave number k incoming from the 
left (right) of the point interaction location. Hereafter, 
subscripts (superscripts) for direction quantum numbers 
indicate that the corresponding r's and t's are those for 
QW (continuous scattering) systems. For the most gen- 
eral zero range potential, we have that (see, e.g., Ref. 
(i?! for a full discussion) 



r 



t 



Cj ± ik{dj 



—Cj + ik{dj + Gj) + bjk^ 

2ik exp[±i9j] 
—Cj + ik{dj + Gj) + bjk^ ' 



6 fc^ 

exp[zLikej], 



(8) 



where a, dj — bj cj ~ 1 , with aj , bj , cj , dj , e j real and Oj S 



[0, 27r) |48|. Equation dH]) satisfies the relations in Eq. © 
and also to 

)(fc) = rf ^*(-fc), if ^(fc) = if ^*(-fc). (9) 



Furthermore, for bj 



0, they become independent 



on k (up to the phases for the r's) and Eq. ([SJ assumes 
the same form as Eq. ([SJ. 

Now, let us set m ~ h ~ 1 (so p — k), define r = 
L/vphasc with Uphaso = p/2, and for convenience take 
L = I. Then, we can make a direct association between 
the quantum walk one step evolution operator U and the 
continuous system propagator Z^(t) , mapping J7|^'(0)) = 

tou{T)\m) ^mr))- 

To concretely establish the correspondence, we start 
with the simplest situation of a fully biased quantum 
walk, i.e., one which always evolves to a same direction. 
We thus assume pj = 1 and 4>j — for any j, from 
Eq. (|4]) leading to J7 = exp[z7] St- Such case presents a 
close parallel with a quantum particle propagating freely 
along the line. In our generalized Kronig-Penney lattice, 
a free particle is trivially obtained by setting all the reflec- 
tion (transmission) amplitudes equal to (1), so that the 
time evolution is U{t) — exp[— i(p^/2)t], with p\p) — p\p) 
for \p) the moment eigenstate. Hence, we have a direct 
mapping between the complete biased quantum walk dy- 
namics and the evolution of a free particle on the line for 
t = T. The equivalent quantities are listed in Table I. 

Next, we consider that in Eq. ([5]) for any j 7^ we have 
Pj — 1 and (pj =0, and for j = we have arbitrary p 
and phases. Also, we assume as the quantum walk initial 
state 



3=0 



*(0)) = ^ exp[zj7]| + l,j) 



(10) 



so that Pj{n = 0) = for j > 0. Then, applying n times 
the evolution operator, Eq. (|4]), to |^'(0)) we get (with 
r = 7-+i,o and t = <+i_o) 
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TABLE I. The correspondence between quantities in the complete biased quantum walk and in the ID free propagation 
Fully biased ID quantum walk Free quantum propagation on the line 

(7 = exp[i7] St W(r) = exp[-j(pV2)T], t = i/wphaso = -f'/(p/2), p = fc, L = l 

l*(0)> = ^ j:',ZtZ expfoTll + 1, j> l'I>(0)> = b) = ^ dxeMip^]\^) 

|*(1)> = (7|*(0)) = |*(0))^ l-l-M) = W(r)|$(0)) = expHpV/2]|$(0)) = expHp]|$(0)) 

J P 

Strictly speaking, exp[i7] in the definition of U in Eq. Q is not necessary and then the equivalence would be complete (here, in each 
step it is up to such global phase). Nevertheless, if we take as an initial state |^'(0)) = instead of states like Eq. l llOII . exp[i7] 

becomes useful to study interference in more general topologies, as in Sec. VI. The important point, therefore, is that we have just a 
phase difference, so not compromising any parallel between the systems time evolutions. 



|*(n))=t/"|*(0)) = --={ J2 exp[ij7]| + l,j)+r ^ exp[-zj7]|-l,.?)+<^exp[ij7]| + l,j)}. (11) 



Now, defining |*scat.) = lim„^+oo !*("-)), one finds 

3=0 ]=-l 



;=! exp[u7]| + l,j) +r ^ exp[-ij7]| - 1, J> + < exp[u7]| + 1, (12) 



l^scat.) ^ , 

2 = -oa J = -oo J=+l 



r 



Note that J7|^'scat.) = l^scat.), so it is a stationary state. are identically null, namely, rj^^ — and t^^^ = 1 (j ^ 0) 

An equivalent situation for the generalized Kronig- ^(+) ^ ^^^^^ ^(+) ^ ^^^y ^^^^ ^^^^ 

Penney lattice is to assume that all r s but one (the re- ^^^-^g ^^j^^^^^ ^ n^c^A^^t from the left reads 

nection amplitude for the point interaction at the origin) 



1 r /"^ /■+°° 1 

|*i^scat.) = —j=\ ] dx exp[ipx]\x) + r{p) / dx exp[—ipx]\x) + t{p) / dx exp[ipx]\x) >. 
V27r I- J-oo J-oo Jo 



(13) 



I 

Comparing Eqs. (fT2|) and ([T3l) . it is evident the corre- a direction to proceed, whereas at the other sites the 

spondence between the two situations. direction is always maintained, with pj — 1 and (fjj ~ 

We can go further, considering that only at two sites for any j ^0, 1. Thus, repeating the same calculations 

the walk can "choose" (from r±ij and i±i,j, j — 0, 1) for the initial state |^'(0)) of Eq. (ITOl) . we get 

I 



j = -oo i=-oc J =+2 

+a| - 1,0) +5exp[i7]| + 1,1)|, (14) 



where 



, , . L t+i,or+i,iexp 2*7 t+i,o , . 

= ^^+1,0 + <-i,oa, i = t+i.i&, a=- -— r, b=- -— -. (15) 

1 - '"+i,i7'-i,o exp[2i7j 1 - r+i,ir_i,o exp[2z7j 

This expression should be compared with that for the associated problem of two general point interactions located at 
X = and a; = 1, whose scattering state (incoming from the left) is given by 

\ r fO rO r+oo 

l^scat.) — \ / dx exp[ipx]\x) + r{p) / dx exp[—ipx]\x) + t{p) / dx exp[ipx]\x) 

V27r ^ J-oo J-oo Ji 

+a{p) / dx exp[—ipx]\x) + b{p) / c?xexp[zpx]|a;)|, (16) 



for the coefficients r(p), t(p), a(p) and 6(p) obtained from Eq. ([T5|) through the substitutions r±ij — >• A'^^p), 
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t±ij — )■ t^j^\p) and "f ^ p. Once more we find a di- 
rect association between the two cases. 

By repeating this procedure of "turning on" more and 
more sites in the quantum walk and zero-range potentials 
in the Kronig-Penney lattice, one realizes that their rela- 
tion is indeed direct. The one-to-one mapping is a simple 
identification of quantities in the two cases. The direc- 
tion coefficients r^ j and i^j" at each site in the quantum 



(a) 



„(±) 



and 



walk corresponds to the scattering amplitudes 

tj^"* of a point interaction in the Kronig-Penney model. 
The quantum number j is associated to the appropriate 
position eigenvalues x/L, likewise for a with respect to 
the signal of p. Lastly, the SQW single step evolution 
U (up to the phase exp[i7]) is akin to U{t — t) for the 
continuous scattering system. 

We finally note that we have discussed the mapping as- 
suming a scattering scenario, with the QW initial state 
given by Eq. (|10p . However, we also could start with 
an initial state localized in some region of the quantum 
walk lattice and an initial wave packet localized in an 
equivalent region of the generalized Kronig-Penney lat- 
tice. Then, by applying the respective time evolution 
operators, again we would find a direct association be- 
tween their dynamics: the multiple scattering of the wave 
packet in the Kronig-Penney lattice resembling the pro- 
liferation of paths (e.g., see the example in Sec. II) in 
the quantum walk. So, the correspondence between the 
two systems is complete and not restricted to the type 
of initial state assumed. This fact becomes more evident 
from the Green function approach next. 
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FIG. 3. (a) If under a particular instance, the QW relevant 
dynamics is restricted to the states \j\ < M; (b) then, effec- 
tively the system can be described by a finite Hilbert space 
"lattice"; (c) whose mapping leads to a finite set of general 
point interactions on the line. 



as illustrated in Fig. 3, in both situations the relevant 
dynamics for the QW is related just to a segment of the 
infinite generalized Kronig-Penney lattice, encompassing 
2M -I- 1 (for M equal to J or N) point interactions. So, 
in such instances effectively one needs to deal only with 
finite lattices. 



B. Quantum walks and finite lattices 

To calculate the exact Green function in the case of an 
infinite generalized Kronig-Penney lattice is a difficult 
task fssl. However, a key aspect in solving QW through 
the proposed mapping is that in almost all situations of 
interest, the original system can be mapped to a finite 
lattice - a limited number of point scatters along the line 
- and not to an infinite comb of zero-range potentials. 

For example, let us assume that the quantum walk ini- 
tial state |*(0)) is, say, either | - 1,-1) or | -|- 1,-|-1), 
thus localized about and leaving from the origin. Now, 
suppose we shall discuss any quantity for times no longer 
than n = N, or for contexts where the dynamics never 
takes the system beyond the sites j — ±J, J > 0. Ex- 
amples are: (a) to determine the probability to be at the 
state j (i.e., to calculate |(j, iT|^'(n)) p) for n up to n = N; 
and (b) to obtain the probability for the walk to reach 
for the very first time a "distance" j = J from the origin 
(j — 0) at times n = 1, 2, . . ., known as the first passage 
time problem in classical random walk theory |49] . 

For (b), any evolution leading to (j,CT |f7" |*(0)) ^ 
> J, arbitrary n) has no interest for the problem 
solution [l^. In (a), after N steps the initial state has 
spread at most a distance |j| — N from the origin. Hence, 



C. The finite lattice Green function and its 
relation to the original quantum walk problem 

Once the quantum walk dynamics one shall study is 
mapped to an appropriate (finite) generalized Kronig- 
Penney lattice, the next step is to calculate the Green 
function for the latter. Based on certain techniques [H^, 
i51.] . the way to do so has been developed in Here 
we just summarize the main steps (for details see [38l|b 

Suppose a particle of energy E = for which 

G{xf,Xi;k) denotes its energy- dependent Green func- 
tion. The initial and final positions, respectively Xi and 
Xf, are arbitrary points along the ID lattice (e.g.. Fig. 
4). Then, the exact G (up to a factor {ik)~^, unnecessary 
for our purposes) reads |38l . [soj 

G{xf,Xi;k) = '^Ws.p.exp[iSs.p.{xf,Xi;k)]. (17) 

s.p. 

The sum is performed over all possible infinite "scattering 
paths" (s.p.) starting and ending at the points Xi and Xf. 
For each s.p., the classical action is written as 5s. p. = 
kLg.p., with Lg.p. the s.p. total length. The pre-factor 
amplitude (or weight) Wg.p. is given by the product of 
the quantum coefficients gained each time the particle is 
scattered off by a given contact potential along the way. 
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FIG. 4. (a) For a finite Kronig-Penney lattice of six gen- 
eral point interactions, and specific end points Xi and Xf, 
a representative "scattering path" composed by stretches, 
(i)-(vi), of straight trajectories. It has a total length of 
Ls.p. = 11 + (-2 - X,) + (1 - Xf), {L = 1). (b) The 7^'s 
and T's in Eq. (|18p are the resulting composite reflection and 
transmission amplitudes for sets of point interaction poten- 
tials, as illustrated. 



To illustrate a typical term in Eq. (|17l) . we consider 
a lattice with six point interactions equally spaced by 
L — I. Taking as the end points —3 < Xi < —2 and 
< Xf < 1, a representative scattering path is depicted 
in Fig. 4. For such s.p., the particle starts at Xi, goes 
to the right, reflects from x ~ —2, moves to the left, 
reflects from x — —3, and then goes to the right, tun- 
neling all the potentials until reflecting from x = 2. In 
this part of the trajectory - stretches (i), (ii), and (iii) 
in Fig. 4 (a) - the partial weight is W(i)+(u)+(ui) = 

ri"^^ tL^2^ iL+i^ t'^o^^ t^"^^ r%\ From x = 2, the particle 
travels to the left, is transmitted through the potentials 
at X = 1 and x = 0, and then is reflected by the point 
interaction at x — —1. Next, it travels to a: = 1 (tun- 
neling the potential at the origin) , suffers another reflec- 
tion, and finally gets to the end point Xf. In this part 
of the trajectory - (iv), (v), and (vi) in Fig. 4 (a) - 
the amplitude is = t'^+i t'f '' r^'^h^'^'' r^:^^\ 

Hence, the total pre-factor weight for this particular s.p. 



is W, 



s.p. 



Wi 



X W( 



{iv) + (v) + (vi) 



The scattering 



1 



path length is simply Lg.p. = 11 + (- 
as readily seen from Fig. 4 (a). 

To obtain (S* in a closed form, one should classify and 
to sum up (c.f., Eq. pT|) ') all the infinitely many possible 
trajectories of the kind exemplified above. Fortunately, 
it always can be done by regrouping the infinite sets of 
trajectories into a finite number of classes [H, [52|. Fur- 
thermore, as proved in fsoj, these classes form geomet- 
ric series, allowing their exact summation. For instance, 
from such procedure the correct Green function for the 
system in Fig. 4 can be calculated, leading to [H, [13] 



G{xf,x^;k) = [1 _ TZiTZ+][l ~ Tl Tl ] - T^T Tl\Tl ("^^Pt"^^^'] + exp[ifca;j]) ( exp[iA:a;/] -I- 7^r exp[-ifca;/]y (18) 

In the above expression, the TVs and T's are effective reflection and transmission amplitudes, resulting from groups 
of zero range potentials as depicted in Fig. 4 (b). They are explicit given by 



T^-i = r^_^ exp[6ifc]. 



t^^l t^^l e xp[4»fc] 



7^, = r^+^ exp[2zfc] + tl^ 

1 - r\_^ r)^2 exp[2ik] 

(rL+^ - (r^iVi+^ - r^+^ exp[2ifc]) 1^:72^1^^] exph2ifc] 



7^+ = rL"t,^exp[-4ifc] 



r+ = 



rfr^V^+^^ exp[2ifc] + (r'fiVLV - ^o^^ exp[4ifc] 

''-1 "-O 



1 - (^rL^a'^'-V + ''-"i^^'o^^) exp[2ifc] -I- (rL"iVL+' - iL^i^iL+i') ^-q^^ exp[4ifc] ' 
TZ- = TZ+ exp[4jfc], where in TZ+ we exchange (+) i — > (— ) and j = —2 i — > j = 0, 

T- =T+, where in 7+ we exchange (-I-) < — > (— ) and j ~ —2 < — > j = 0. (19) 



r 



Here, two points should be emphasized: (a) the energy domain G is given by a sum over all the possible multiple 



scattering processes suffered by the particle; (b) each s.p. 
in the series Eq. (|17p represents a trajectory in which 
the particle spends a time of roughly t ^ nr, for n the 
number of scattering along the path (e.g., n = 12 in the 
example of Fig. 4). 

In the study of QW, common questions are related to 
the system state, say, after evolving n steps. But from 
the (a)-(b) above, such information is fully contained in 
the series representation of G. Therefore, since the corre- 
spondence between QW and generalized Kronig-Penney 
lattices is straightforward, we can readily associate each 
term in Eq. ()17p to a possible evolution of a quantum 
walk (e.g., that in Eq. ([T])). Moreover, such terms can 
be viewed as a Fourier decomposition of G. Given that 
the Green function Fourier transform is the time domain 
propagator, an individual term in Eq. ()17|) . when prop- 
erly mapped, represents then a possible path for t — n 
time steps in the quantum walk. 

Finally, depending on specific QW quantities we shall 
calculate, in practice a simple inspection and selection of 
paths in the expansion for G will suffice. However, for 
larger n's and more complicated topologies (Sec. IV), 
it may be cumbersome to deal with individual terms in 
Eq. ([TT)) . Fortunately, one can make the Green function 
a systematic protocol for QW by introducing the path 
and step operators. As we discuss in Sec. V, they are 
useful tools to extract any information about the system 
directly from an already summed closed expression for 
G. 



IV. EXTENSION TO ARBITRARY 
TOPOLOGIES 



(a) 



1 1 



r' 'Jx ^'■'yj^_< ; I ■' "^x ''-ly'/ I ' '^x'-'y'/j; I ' "^x'-'y'/ 

|2, a -uTfyorij > |2, (j,, jy)>! l3:a,, jy)> 




FIG. 5. Examples of graph structures, which generalize ID 
QW. For SQW, the states (as illustrated) are defined on the 
bonds joining the different sites j. (a) All the sites have a 
same number of first neighbors in a regular topology, (b) For 
an irregular structure, such number depends on j. 



QW can be defined in arbitrary topologies [15|, i.e., 
for general graph structures [TJ. The scattering formula- 
tion is then obtained through a direct extension of the 
construction in Sec. II [1, [l,13l- 

Suppose an undirected simple network ^] of sites con- 
nected by bonds (examples in Fig. 5). Its topology repre- 
sents the particular Hilbert space arrangement in which 
the quantum walk dynamics takes place. Like the ID 
lattice, each bond joining two neighbor sites - say j and 
j' - is associated only two basis states, one incoming to 
J and other to j' . For instance, for the bond joining 
{jx - l,jy) to {jx,jv) in Fig. 5 (a), we have |1, {jxjy)) 
and \3,{jx — l,jy)), whereas for the bond connecting j to 
f (j to j") in Fig. 5 (b), we have |3,j) and {\2,j) 
and |2, j")). But contrary to the ID case, the possible 
"directions" (bonds) to get to a site j from its first neigh- 
bors may depend on the specific j. Thus, the quantum 
number Cj assumes the values 1,2,..., Kj, with Kj the 
coordination number (valence) of site j. In more regular 
structures Kj — K regardless of j (e.g., K = A\n Fig. 5 
(a)). 

The construction of the time evolution operators is 
discussed, e.g., in Refs. [1, Here we just outhne 

the main ideas following Ref. Q- First, one needs to 



characterize the lattice topology, namely, to specify for 
any j the sets: (a) Sj = • . . , j^^^^} of ah the 

Kj sites which are first neighbors of j (e.g., in Fig. 5 
(b), S, = {j', /"}); (b) AG = , a,(2, , . . . , a } 
for aj{n) the direction quantum number for the state in- 
coming to site j*^") through the bond joining j^^"^ and 
j (in Fig. 5 (b), Nj = {aj^ = l.cjjn = 2,aj,„ = 1}); 



and (c) Bj = {(i (i),(t 



,a-.(Kj)} with fT („) the di- 



rection quantum number for the state |(Tj(„),j) incom- 
ing to j along the bond connecting j and j*-"-* (in Fig. 



5 (b), Bj = {a, 



3 : bond j-j' , (T,(2) 



bond j-j", djm = 1 : bond j-j'"}). 

Second, one defines the one step time evolution U in 
terms of its action over each basis state |tTj, j), or (with 
(Tj ~ dj{i) G Bj and UjH) the corresponding element in 

AG) 

U{l)\crj,j) = exp[i7](^r^^.^^.j|crj-(,),i**^) 

n=Kj 

+ E i-.^,(.,,.k.<">'J^"^))- (20) 
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Finally, the r's and t's are chosen such that for any j Sj\k) = Sj{~k)), yield (where i,l,n = 1,2,..., K^^ 
the Kj X Kj matrix Tj (of elements [r^Jo-g- — ra-a-.j and 



[^j]cr'a — tc'a.j, toT both <T 7^ (t' in {1,2, Kj}) is 
unitary. This makes U also unitary [7], establishing SQW 
in arbitrary topologies. 



The usual, i.e., continuous in time and space, quan- 
tum mechanical dynamics on network structures (known 
as quantum graph systems [13]) is likewise a generaliza- 
tion of the evolution in a ID lattice with zero range po- 
tentials [H, [11] . It is obtained by properly matching the 
solutions of a ID free [111 Schrodinger equation in the dif- 
ferent arms (bonds) at the vertices (sites). For this end, 
one assumes for each j a matrix Sj{k) (see below), which 
describes how an incoming plane wave of wave number k 
is scattered off at the vertex j. So, any j can be viewed as 
a general point interaction, but connecting i^^-directions 
instead of the common two (left and right) on the line. 
Furthermore, if for all j, SjSj'' = Sj^Sj = 1, the result- 
ing dynamics is unitary, conserving flux probability. 

Quantum graphs can be solved through a Green func- 
tion approach ^57,]. In fact, it has been shown 39] that 
the exact G is also given by Eq. ((T7]), where now the 
scattering paths are all the possible trajectories along 
the network, starting and ending at the points Xi and 
Xf (located in arbitrary arms of the graph). The W^'s 
are the quantum amplitudes gained along the s.p.'s due 
to the scattering at the different sites. Finally, the men- 
tioned procedure of classifying and summing up different 
classes of s.p.'s still holds in this case IBS]. So, we always 
can write the exact G as a closed analytical expression. 

Summarizing, QW in general networks are direct ex- 
tensions of QW in ID exactly in the same way than quan- 
tum graphs are natural extensions of Kronig-Peney lat- 
tices. Therefore, it is easy to realize that our previous 
mapping between the two types of systems in ID remains 
valid in arbitrary topologies too. 

Lastly, to define the reflection and transmission scat- 
tering amplitudes in a quantum graph - and to associate 
them to QW coefhcients - we consider the same label- 
ing used to characterize the lattices topologies. Thus, for 
Bj — {(Tj(i) , (Tj(2) , . . . , ^jiKj ) }, the matrix elements of Sj 
are 



3 y^j — ^] 



[Sj]ii 



(5. 



)) 



In Eq. dH 



[Sj]ii 



(21) 

') is the reflection 
(transmission) coefficient for the particle incoming to site 
j from bond o'j(i) to be reflected (transmitted) to bond 
CTj(i) ((Tjco). The unitarity of the S*-,- 's plus the symmetries 
of the Schrodinger equation for real potentials [SSi] (i.e.. 



tf^<*'''^<"^(fc)=tf^<"*^«^*(_fc), 



(i-)) (Sj(i)Sj(i)) 



j j J 



(22) 



The above are natural generalizations [3^,[53 of the usual 
relations for the scattering coefficients (c.f., Eqs. and 
^) of a point scatterer on the line. Note also that if we 

impose time-reversal mvariance, i^- ~ j 

Hence, the direction coefhcients in a quantum walk, the 
Fj's, are in one-to-one correspondence with the scattering 
matrices S'j 's in a quantum graph system. 



EXTRACTING INFORMATION FROM G: 
THE STEP AND PATH OPERATORS 



From the previous results, it turns out that the exact 
Green function, Eq. (jl7p , is actually the generating func- 
tion of all the possible walks leaving from and arriving 
at the bonds corresponding to Xi and Xf, respectively. 
So, any individual or group of QW paths are obtained 
through proper manipulations of G. 

In this way, more simple tasks like to determine all 
the trajectories for |^(0)) — \a,j) evolving, say, only two 
times steps (n = 2), can be done by identifying particular 
terms directly in the G series representation, Eq. (|17p . 
However, the huge proliferation of paths in instances such 
as to find certain trajectories connecting two bonds very 
far apart, or resulting from high values of n, makes the 
full series expansion difficult to deal with. In such cases, 
a better approach is first to sum up the series [s^ (using 
the already mentioned procedures in the literature [3^, 
m, [sOjUII to get expressions like Eq. ([18])) and then to 
extract the sought information from G with the help of 
the two operators described below. 

The first is Sn , yielding all the paths of exactly n time 
steps. To define Sn, note that any walk state gains a 
factor exp[i7] at each time step (see Eq. ([^0]) ). From 
the mapping, such factor is equivalent to z = exp[ifcL] 
in the (continuous) quantum graph problem. So, let 
us set Gz as G in Eq. (fTT]) . but with the substitution 
exp[?fcL] — z, and for which the scattering amplitudes 
are identified with the appropriate quantum walk coeffi- 
cients Fj 's (Section IV). Thus, one finds that if the step 
operator Sn acting on Gz has the form 



Sn = 



1 d" 



n\ dz" 



(23) 



z=0 
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I " 1 2 . 

then SnGz gives the total probabihty for the quantum 
walk to leave the bond Xi and to get to the bond in ex- 
actly n steps. We should mention that such construction 
has already been proposed in 0, 13] to treat scattering 
problems. Considering the Green function approach, we 
see that S'„ can be applied in more general contexts. 

The second is P-p, which extracts from G all the paths 
with specific trajectory stretches V . Any quantum walk 
s.p. can be described by the sequence of coefficients r 
and t it gains along the way (cf., Eq. ([7])). For instance, 
consider Fig. 5 (b) and assume n = 6 applications of U 
to the system initially at One possible sequence 

of successively visited states during the evolution, thus 
representing a possible path, is: 



|1,jH|1,/)^|3,j)^|1,/)^|3,j-)^|2,/')^|2,j). 



Here W = hij riij> r^aj riij> t32,j r22.,j" is its prob- 
ability amplitude, which can be rewrite as W = 
(^33j)^ ('"22j")^ (^iij)^ (^32j)^- Thus, any tra- 
jectory (or part of it) can be represented hy V — 

{(ai,naJ,(a2,naJ,...;(/3i,n/3i),(/32,n/32).---}, i-e-, by 
the set of indexes and exponents associated to the quan- 
tum coefhcients r and t of the path stretch (with a and 
/3 standing for a'a",j). In our example, ai = 
"2 = 3 3, j; as = 22, j; /3i = 1 (32 = 3 2, j; = 2; 
= Has = np^ = Ufj^ 1- 

Now, by properly choosing Xi and x f (which obviously 
depends on V, see Sec. VI) we get that G-p = P-pG is 
a sum ~ in the form of Eq. p7p - but containing only 
paths whose parts of their trajectories are given by V. 
The explicit expression for P-p is 



(24) 

which must act on the Green function as the following: 
first one performs all the indicated derivatives; second, 
one sets the coefficients and equal to zero; finally 
one multiplies the resulting expression by the coefficients 
(ra)"<»'s and (t/3)"^'s. 

If we shall select just a path which is itself entirely 
represented by V , then in the above definition we simply 
change |r„=o a-nd |t^=o by \r=o and \t=o, with r and t all 
the quantum amplitudes in G. 

Lastly, we note that for an initial state being the su- 
perposition of N basis states, |^'(0)) = X^o-j Cct.jIc, j), we 
must consider N Green functions, each with a Xi corre- 
sponding to the bond of |cr, j). So, in any calculation, the 
contribution of each of these G's should be weighted by 
the associated factor Caj. 



VI. AN EXAMPLE: A DIAMOND-SHAPED 
GRAPH 



Finally, to illustrate some features of our Green func- 
tion approach, we discuss a quantum walk for the topol- 
ogy depicted in Fig. 6. We assume complete general 




FIG. 6. (a) A graph composed by a diamond-shaped region 
(sites A, B, C, D) attached to semi- infinite leads (sites j < — 1 
and j > 0, for which rj — and tj = 1). (b) The a labelhng 
for A, B, C and D. (c) For n — 5 time steps, schematics of all 
possible s.p.'s first entering into the diamond region through 
the superior bond (those entering through the inferior bond 
are simple specular images of (c)). 

coefficients (observing Eq. ([22])), in the diamond region 
- sites A, B, G and D - and free evolution, rj — and 
tj = 1, in the leads region - sites j < —1 and j > 0. We 
mention that this system, in the case of r^i = ro = —1/3, 
tA = to = 2/3, rB ~ rc = and tB ~ tc ^ 1, has been 
studied in Ref. . Such particular values represent the 
so called Grover coins (see, e.g., [3]). 

Consider this quantum graph for the initial and the 
end positions Xi and Xf, respectively, in the bonds i and 
/, Fig. 6 (a). The exact Green function reads 

G{xf,Xi;k) = Texp[ik{xf — Xi)], (25) 

with T the global transmission amplitude resulting from 
the multiple s.p.'s which cross the diamond-shaped re- 
gion. By using the procedures in i39j , one gets 

P+/g {P-/g) represents the contribution of all the s.p.'s 
which initially enter the diamond region through the su- 
perior (inferior) arm. They are given by (ct = ± and 
a = —a) 
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5 = { 1 - (r+ - t- +,A + i?+ + r+ +^a) exp[2z7] } { 1 - (r_ + t+ + i?_ _ r_ _^^) exp[2z7] j 

T+o = t+^^B 1^+0, _D +?'--,c {t+-.Dt-[),D ~r^^,Dt+o,D^ exp[2i7] | exp[i7]//, 
T+- =t+ _,_B t+ -,D t- +,c exp[2i7]/ /, 

R+- = r++,B + <+-,B|t-+,B?'++,_D + t-+,Br- -,c (t-+,Dt+-,D - r^+^oT- -,D^ exp[2«7] | exp[2i7]//, 
/ = ( 1 - J"- -.5 +,D exp[2i7] ) ( 1 - r_ _,c r_ exp[2i7] ) - <+ t- +.D r- -.b -,c exp[4i7] 



r_o = T+Oi T-+ =T+-., R = R+ + , where in all the r.h.s. terms we raust exchange B i — > C. (27) 



For Eqs. (P5 1) - (P7)) we have already used the mapping, 
writing them in terms of quantum walk quantities. 

The amplitude T simplifies considerably if for any cr, 
a', we have tfjfj' x — tx a-nd ra-a-,x = with X = 

I 



I 

A, B,C, D. Furthermore, if the coefficients for the sites 
A and D and for the sites B and C are set equal, namely, 
TA = ro, tA = to, TB = rc and Ib = tc, Eq. dUIl-dlfl) 
yields 



r 



2*4 tB 



l-2{tA + ta) TB exp[2i7] - {Ia + taY it 



;) exp[4i7 



exp[3i7] 



(28) 



For the particular Grover coin values, we get from Eq. 
([28|) T = 8exp[3i7]/(9 — exp[4i7]), in agreement with 
Ref. [44] as it should be. 

We emphasizes that T is given by a sum over all the 
possible paths leaving i, going into the diamond region, 
and finally leaving to the bond /. So, the probability for 
i — / in exactly n steps can be obtained by applying 
the step operator to 7^ = T|cxp[i7]->z- Such type of cal- 

I 



^5 = *0+,a| 

+ta-,A{ 



I 

culation is useful because it gives the hitting time 
[Slj . i.e., the probability for the walk to reach a given 
state \iy,j) from |cr',j') as function of n. The present 
Green function approach allows to obtain hitting times 
in a rather direct way. To exemplify this, we first con- 
sider the most general case, Eqs. (I^ - (P7)) . and select all 
the paths reaching the bond / in five time steps. Then, 
hh = S^Tz, reads 



t+-^Br++^Dr--^B + r++^Br++,At+-.B i+o,D + 



-.B < 



,D r--,c + r+ +.B t+ -^A t+ -^c 



t+ -,C t- +.D r- -.B + r+ +,c t- +,A t+ -,B 



o,d| 
(29) 



r 



which represents the eight possible trajectories with n 
5, Fig. 6 (c). 



Certainly, in more symmetric situations analytical re- 
sults are easier to obtain. Indeed, for the case of Eq. 

I 



28| . Tz can be casted as 



r. = 



tB +rB , tB - TB 



z+ 



(30) 

where z± = (iis - rs)/[(tA + 7'A)(t| -r%)]. Hence, for 
\h„\^ = !(„!)-! (5"rj9^")z=oP, we find that 



\K\' = \t\{tA + rA)^'^-'^/'~'\'>c 



\{tB + rB)("-i^/' - (-l)("-i)/'(iB - rB)("-i)/2|5 



if n is odd, 
if n is even, 



(31) 
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is the probability to cross the diamond region in exactly n 
steps. |ft-„p = for n < 3, since at least three time steps 
are necessary to leave the bond i and to arrive at bond 
/. Also, from a direct inspection in Fig. 6 one realizes 
that it is not possible a transmission for an even number 
of steps, a result explicit in Eq. ([3T|) . li ts = tc = 0, 
obviously the system never gets to the right lead and 
Eqs. (pS)) and (I5T]) are identically null. Finally, if we 
consider tb — and Ib = 1, we find P„ 7^ only for 
n = 3 (mod 4). Furthermore, assuming tj^ = 2/3 and 



TA = —1/3, we recover the result in Ref. [13], namely, 
|/i„|2 = (8/9("+i)/4)2 for n = 3 (mod 4) and \K\^ = 
otherwise. 

As discussed in the previous section, specific T^'s are 
obtained from the Green function by means of the path 
operator P-p- For instance, suppose we shall select the 
trajectories directly crossing the diamond region, i.e., 
transmissions through A, then through B or C, and fi- 
nally through D, with no multiple reflections. In this 
special case, the path operator is 



Pv = 



d 



1! dt+-.B 



t=0{t^to 



+ c+ a,D) 



1! dt 



(32) 



t = 0(t#to-,A,t-0,£3),'-=0 



leading to 

P-pT ^ {tQ+,At+-.Bt+Q,D+tQ-.At+-.ct-o,D) exp[3i7] 

(33) 

Note that Eq. I\'3'2\i is in a simpler form than the general 
definition, Eq. ([24|) . This is so because we have used 
the fact that Pp acts on a transmission Green function. 
Indeed, there is no need to perform derivatives as, e.g., 
toa-.A (d / dtci a-,A)\ta „_A=o- Thus, the path operator is con- 
I 



siderably simplified if we choose suitable configurations 
to calculate G. 

We can think of more general paths, namely, to cross 
the diamond region in a total of n = n+ + n_ steps, but 
for exactly n+ (ri_) steps taking in the superior (inferior) 
arm, i.e., in the bonds A-B and B-D {A-C and C- 
D). If for simplicity we assume that for each site X 
{X = B or C), all the t's and r's are equal, regardless 
the quantum numbers c's (as in Eq. ([M])), then the 
mentioned operator, to be applied to T, is written as 



Pv 



E 

(1) I 1 (6) 1 (1) I 



(1) (2) (3) (4) (5) (a) 

{t+.,Ar- [t^-^aT- {tcT- [tcT- {t+-,Dr- (r_„,z,)"- 



,(2), 



,(3), 



,(4), 



,(5), 



,(6), 



/ 9"+ \ / a"+ \ / a'"+ \ / c'"+ \ / a"+ \ / a"+ \ 

V „ ) [ „ n!f' ) [TTFJ \~^) V „ n<:'> ) [ „ ni^' ) 



(5) 



I dt 

5"- 



dr 



(6) 



( „ n'^' ) ( „ n<=' ) ( „ n'^' ) ( „ n*^' ) ( „ ) ( „ n'«' ) 



(9< 



dr 



dt 



dt: 



dr 



-.D 



t=0(t#to±,A, t±o.n). r=0 



(34) 



Although the above expression may seem rather cum- 
bersome, it is amenable to work with by using algebraic 
manipulation softwares (what we have tested by investi- 
gating different situations; results will be reported else- 
where). In certain instances, nevertheless, the calcula- 



tions can be carried on straightforwardly. For instance, 
consider all the paths which get to right lead only pass- 
ing through the superior arm. Furthermore, assume that 
among them, we shall select those tunneling the site B 
exactly n times. In this case, the path operator is simply 



Pv = 



{tBT d" 



dti 



(35) 



t=0 (t 7^ to t+o,D), ''=0 (r 5^ r+ +,A, r-++,D) 
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For n even its results in P-pT = 0, and for n odd in 



p^j- ^ ^o+.A (tB)"(r++,D)"^ t+ 0,1) exp[(2n+ 1)^7] ^.gg^j 

(1 - r++,A?'sexp[2i7])-T^(l - rsr++, ^ exp[2j7])-T^ 



VII. CONCLUSION 



By means of appropriate mappings to systems for 
which the exact G can be calculated, quantum graphs, 
we have obtained closed and general expressions for SQW 
Green functions. Furthermore, the procedure allows to 
discuss complete arbitrary topologies and position depen- 
dent quantum amplitudes (60i] . 

By introducing two operators, namely, step and path 
operators, we have shown how to extract from G any 
relevant dynamical information about the system. For 
instance, one can exploit particular paths in a quantum 
walk as well as to obtain the contribution of orbits of spe- 
cific time length n. Such possibilities have been exempli- 
fied in details for a diamond-shaped graph structure. 

An important fact, not explored in this contribution, is 
that our formulation naturally allows the introduction of 



energy (or equivalently, wave number k) dependent tran- 
sition amplitudes. In QW context, such k could be faced 
as an extra inner variable. Since different walks may have 
different /c's, using the Green function approach, then 
one could address the case of collective QW. A complete 
study of energy dependent SQW will be the subject of a 
future work. 

Finally, we have discussed G only for QW scattering 
formulation. Nevertheless, as already mentioned in the 
Introduction, the SQW and CQW are unitary equivalent. 
Moreover, CTQW are associated to CQW. So, the Green 
function framework for SQW can be extended to such 
other constructions as well. 
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